Eccentricity and elliptic flow in proton— proton collisions from parton evolution 
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It has been argued that high-multiplicity proton-proton collisions at the LHC may exhibit col- 
lective phenomena usually studied in the context of heavy-ion collisions, such as elliptic flow. We 
study this issue using DIPSY — a Monte Carlo event generator based on the QCD dipole model. 
We calculate the eccentricity of the transverse area defined by the spatial distribution of produced 
gluons. The resulting elliptic flow is estimated to be about 6%, comparable to the value in nucleus- 
nucleus collisions at RHIC and the LHC. Experimentally, elliptic flow is inferred from the azimuthal 
correlation between hadrons, which receives contributions from collective flow, and from various 
other effects referred to as "nonflow". We discuss how to identify in experiments the signal of flow 
in the presence of large nonflow effects. 
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I. INTRODUCTION 

Elliptic flow is one of the most important phenomena 
observed in ultrarelativistic nucleus-nucleus collisions 0- 
3]. An Au-Au collision at RHIC produces several thou- 
sands of particles. If interactions among these particles 
are strong enough they expand collectively like a fluid, 
and elliptic flow is a probe of this collective behavior Q ■ 
The fluid picture is a macroscopic one, which is gener- 
ally valid for a large system. For a system as small as a 
nucleus, it is an idealization which must be amended in 
order to quantitatively understand experimental data |5| . 
The system formed in proton-proton (pp) collisions is 
even smaller. Yet the possibility has been raised that el- 
liptic flow may be seen in pp collisions at the LHC [fl-flijl . 
This is actually a quite nontrivial problem which can only 
be addressed with a proper understanding of the proton 
wavefunction at high energy_from QCD, whereas most of 
the preceding works [6, 9- 13] are based on rather primi- 
tive models of the proton. [See, however, [7l[l3|.] In this 
letter, we study this issue using a full Monte Carlo (MC) 
model of the collision which implements the BFKL-type 
evolution of structure functions, multiple collisions, the 
partonic shower and the subsequent hadronization. This 
model is briefly described in Sec. ILT1 

An obvious obstacle to develop collective phenomena 
in pp collisions is the low multiplicity of hadrons in the 
final state. This may be overcome by triggering on high- 
multiplicity events. Indeed, it has already been observed 
in the 7 TeV run at the LHC [H,[ll| that the multiplicity 
distribution has a broad tail reaching out to ^ ch > 30, 
and this will be further pronounced in future runs at 
14 TeV. Such high-multiplicity events originate from up- 
ward fluctuations in the gluon multiplicity inside the pro- 



ton and the subsequent multiple gluon-gluon scatterings. 
The fluctuations in the distribution of gluons in the trans- 
verse plane then generates nonzero eccentricity of the in- 
teraction region even in collisions at vanishing impact pa- 
rameter. 1 Assuming hydrodynamic evolution for these 
high-multiplicity events, we estimate the magnitude of 
the resulting elliptic flow in Sec. IIIII 

Experimentally, elliptic flow is not measured directly, 
but inferred from azimuthal correlations between the pro- 
duced particles. These correlations are partly due to el- 
liptic flow, partly due to other effects referred to as "non- 
flow" [l7j. Nonflow correlations are sizable for peripheral 
nucleus-nucleus collisions at RHIC [l8[ , and one expects 
them to be even larger in pp collisions, making it a chal- 
lenging task to disentangle the flow contribution. In Sec- 
tion IV, we stress the necessity to look at higher-order 
cumulants of azimuthal correlations, and suggest how to 
identify flow in light of the experimental and MC results. 



II. THE MODEL 

Our calculations are based on the MC implementa- 
tion of the dipole model developed in Lund [19h24| . The 
dipole model by Mueller [2^, [2(| realizes the leading- 
order BFKL evolution of gluons in the transverse coor- 
dinate space, which is ideally suited for the computation 
of the eccentricity. It is known that the BFKL evolution 
generates large event -by-event fluctuations in the gluon 
multiplicity [27| as well as characteristic spatial distri- 
butions and correlations in the transverse plane (28l . l29l | . 
Both of these effects are important in properly estimating 
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Previous works considered the fluctuation of 'hot spots' I lull and 
'flux tubes' lllll . There the transverse distribution of these ob- 
jects was assumed to be random. In our case the transverse 
distribution of gluons is governed by the QCD evolution. 
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the eccentricity. 

For phenomenology, the original leading-order formu- 
lation is impractical, and over the years there have been 
many improvements of the model. These include the run- 
ning of the coupling, energy-momentum conservation, 
saturation effects, and confinement effects at large dipole 
separations. [For details, see 1!) 22|.| The parameters of 
the model are determined so that they reproduce a few 
observables (such as the total pp cross section) at some 
energy. Predictions can then be made for various other 
observables at different energies without any further tun- 
ing of the parameters. 

The MC code that we actually use in the following, 
called DIPSY (23l . |24| , is the most advanced version by 
the Lund dipole team which has access to all the ex- 
clusive final states. In this framework, a typical high- 
multiplicity event looks as follows: Before the scatter- 
ing each proton develops a cascade of gluons (or equiv- 
alently, dipoles) spread in rapidity and the transverse 
plane. These gluons, mostly soft ones, then undergo mul- 
tiple scatterings. The evaluation of the non-diffractive 
scattering amplitude for the two cascades reduces to that 
for individual pairs of dipoles, allowing DIPSY to decide 
on an event-by-event basis which dipoles interact. It is 
then possible to trace the interacting parton chains back 
from the interactions, and the initial state radiation can 
be identified. All emissions not connected to the interact- 
ing chains are reabsorbed as virtual fluctuations. The ini- 
tial state radiation is then passed to ARIADNE [3(3] that 
further splits the dipoles with timelike emissions. After 
that the dipoles hadronize through the string fragmen- 
tation model in PYTHIA (3ll. l32l|. giving the observable 
final states. 



III. ECCENTRICITY AND THE ELLIPTIC 
FLOW 

In a nucleus-nucleus collision, the participant eccen- 
tricity e par t is defined from the positions (x, y) of partic- 
ipant nucleons within the nucleus [3j: 
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where (• • • ) denotes averaging over events in a given cen- 
trality bin. Hydrodynamic evolution linearly relates e{n} 
and the corresponding elliptic flow i^ jn j measured from 
the n-particle azimuthal correlation [33J. [See, however, 
a recent study [34] which suggests a possible mixing of 
different harmonics due to fluctuations.] An empirical 
formula which works at RHIC is Q 
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and a similar relation between i>2{4} and e{4}. In ([5]). 
(^2/e)hydro ~ 0.2 is the ideal hydrodynamics result and 
the parameter X/Kq = 5.8 fm -2 measures the degree of 
incomplete equilibration. S is the area of the overlap 
region calculated as 
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1.5 is the total hadron rapidity distribu- 
tion. [We neglect the small difference between the rapid- 
ity and the pseudorapidity.] 

In pp collisions, we employ the same formulae (dJ) — (J3J) , 
with the averaging in |2j performed for the "liberated" 
gluons, i.e., those in the initial state radiation. In doing 
so, we apply a rapidity cut such that only gluons which 
are separated from the beam directions by more than 2 
units of rapidity are included in the averaging ([2]). 2 The 
validity of the use of ([5]) in pp collisions at the LHC is 
a priori not clear and requires an explanation. We first 
note that hydrodynamic simulations show that i>2 /e as a 
function of ^ falls essentially on the same curve both 
at RHIC and at the LHC in spite of the difference in en- 
ergy by a factor of about 14 [6|, and both for Au-Au and 
Cu-Cu at RHIC although they differ in size by a factor 
of two [The parameters (t>2/e)hydro and X/Kq in prin- 
ciple depend on the temperature, but their changes are 
small due to the softness of the QCD equation of state.] 
We then recall the general argument that the applicabil- 
ity of hydrodynamics is controlled by the dimensionless 
parameter a = d ^ drj rather than the system size. In 
our simulations, the necessary condition a < 1 is well 
satisfied in a broad N c h range even after allowing for 
some uncertainty in the parameter X/Kq. On the other 
hand, it is hard to imagine hydrodynamic behaviors in 
systems smaller than S ~ I fm 2 which roughly sets the 
border between the hadronic and nuclear scales. We thus 



and the brackets {• • • } denote averaging over the partic- 
ipants in a given event. We shall be interested in the 
quantities e{2} and e{4} defined by 
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2 We could have included only gluons which are centrally pro- 
duced, say, within \r]\ < 1. However, we find it more reasonable 
to average over a wider range in rapidity in order to make up 
for the lack of the final state radiation in the eccentricity calcu- 
lation. Noncentral gluons are connected to the central ones by 
color stings, and the early stage final state radiation from the 
noncentral gluons will contribute to the central energy density 
(and vice versa) before hydrodynamics starts to operate. 



3 



expect that §5$ can be marginally applied for S > 1 fm 2 , 
and this condition is better satisfied in high-multiplicity 
events (see below). 

For the actual evaluation of v 2 , we propose the follow- 
ing slight improvement of ([5]) 
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The reason is that in pp collisions the eccentricity e and 
the area S fluctuate widely even at fixed dN/drj. Equa- 
tions ([7|) and fJ3J) nicely captures this event-by-event cor- 
relation between e and S. Note that it is the squared 
value (^2{2}) 2 (and also (t>2{4}) 4 ) that directly comes 
out of the experimental measurement of flow via multi- 
particle correlations pj|] 
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where fa is the azimuthal angle of the i-th outgoing par- 
ticle and averaging over all pairs (and 4-plets) satisfying 
some cut requirements is implied. In the case of nucleus- 
nucleus collisions, the fluctuations are rather small so 
that ([7]) essentially reduces to the previous formula ([5]). 

We have generated pp events at yfs —7 TeV and 14 TeV 
with randomly chosen impact parameter b, and classified 
events in bins of the charged particle multiplicity N c h- 
The averaging (• • • ) has been taken in each bin. Unlike in 
nucleus-nucleus collisions, in pp collisions the impact pa- 
rameter is not measurable, and there is no simple scaling 
between the centrality and the multiplicity (not even be- 
tween the centrality and the effective area S) because of 
the fluctuations. Still, the majority of high-multiplicity 
events in our MC simulations comes from collisions with 
6«0. 

In Fig. [TJ we plot the results for e{2}, e{4} and (S) as 
a function of N c h within the ALICE acceptance \rj\ < 0.9 
(central detector) at 7 TeV. Events with N c h = 60 typ- 
ically have 12 dipole-dipole (gluon-gluon) subcollisions. 
We see that the eccentricity is 20-40% in the highest 
multiplicity region, similar to the value in semi-central 
nucleus-nucleus collisions. 

If the rescatterings among liberated partons (not in- 
cluded in DIPSY) are strong enough, then the nonzero 
e will give rise to v 2 according to the modified formulae 
0, ©. We plot the results for v 2 {2} and u 2 {4} in Fig.d 
We see that v 2 is more or less constant and is about 4- 
6%. This is comparable to the value in nucleus-nucleus 



collisions at RHIC and at the LHC. We discuss the impli- 
cations of these results below. The Monte Carlo shows no 
significant difference at 14 TeV (results not shown) even 
though high- multiplicity events are then more frequent. 
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FIG. 1: (Color online) Predictions for the eccentricity and the 
interaction area S [fm 2 ] (xO.l) at 7 TeV versus the charged 
multiplicity in the interval |r/| < 0.9 at yfs = 7 TeV. 

It is worth noting that DIPSY predicts that the con- 
ventional definition of the eccentricity 



(ii) 



typically takes a negative value, if the impact parameter 
is chosen in the s-direction. This is in stark contrast 
to the nucleus-nucleus case where the interaction region 
is roughly the geometrical overlap of two colliding nuclei 
so that e s > always. While e s is unmeasurable in pp 
collisions, this still illustrates the fact that the origin of 
the eccentricity is very different from that in the nucleus- 
nucleus case. 



IV. DISCUSSION 

The rather large value of elliptic flow v 2 ~ 6% that we 
have obtained in the flow scenario may seem promising 
at first. However, in practice the observed i>2{2} and 
D2{4} differ from the genuine v 2 by the so-called nonflow 
contribution 



(v 2 {n}) n = v n 2 + S n , 



(12) 



where 8 n is the n-particle correlations not associated 
with flow, such as resonance decays and the back-to- 
back correlations from hard and semi-hard scatterings. 
In nucleus-nucleus collisions, they are relatively innocu- 
ous because they scale with the multiplicity as 
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FIG. 2: (Color online) Predictions for «2{2} and t>2{4} calcu- 
lated from Eqs. (0 and @ versus the charged multiplicity in 
the interval \rj\ < 0.9 at = 7 TeV. The error bars are from 
statistics only. 

In pp collisions, we expect that the nonfiow contribution 
is less suppressed than (|13[) due to various initial and 
final state effects. Most importantly, high-multiplicity 
events often contain several jets and mini-jets, and par- 
ticles within a pair of recoiling jets typically give a large 
contribution cos2(0j — <fij) « 1 to the average There 
are also correlations from the initial state partonic evo- 
lution which is initiated by only a few partons. These 
correlations tend to enhance the nonfiow contribution 5 n , 
making the isolation of t>2 difficult for small values of n. 

Indeed, the ALICE collaboration has found that 
i>2{2,4} decrease slowly with N c h (35[. This slow de- 
crease is also observed in Monte Carlo simulations, mean- 
ing that the scaling (|13|) does not hold for pp collisions. In 



the highest multiplicity events, V2{2} s=* 0.13 HH which 
is twice as large as the flow contribution v 2 ~ 0.06. This 
implies that the two-particle correlation is dominated by 
nonfiow effects. 

This situation in pp collisions necessitates us to look 
at higher order cumulants V2{n} with n > 4 which 
are by definition insensitive to two-particle nonfiow 
correlations. In this regard it is very interesting to notice 
that the ALICE collaboration [35fl also reported the 
measurement of i>2{4} in pp, and found that the right- 
hand-side of (fTU|) is negative. [(«2{4}) 4 ~ — 10~ 4 in the 
highest multiplicity bins.] The same phenomenon can be 
seen in MC simulations like PYTHIA (as done in j35|) 
and also in DIPSY without assuming flow. On the other 
hand, as shown in Figs. Q] and [2J (e{4}) 4 and (^2{4}) 4 
are positive in the flow scenario at large Nch- Leaving 
the origin of negative four-particle nonfiow correlations 
for future work, we conclude with a prediction that, if 
there is flow in the large N c h region, then the fourth 
order cumulant (fTUl) . which is negative in the absence of 
flow, will eventually turn positive. 
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